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We present a microscopic derivation of the effect of current flow on a system near a superconductor- 
metal quantum critical point. The model studied is a 2d itinerant electron system where the electrons 
interact via an attractive interaction and are coupled to an underlying normal metal substrate which 
provides a source of dissipation, and also provides a source of inelastic scattering that allows a 
nonequilibrium steady state to reach. A nonequilibrium Keldysh action for the superconducting 
fluctuations on the normal side is derived. Current flow, besides its minimal coupling to the order 
parameter is found to give rise to two new effects. One is a source of noise that acts as an effective 
temperature Te// = bEvfTsc where E is the external electric field, vf the Fermi velocity, and Tsc 
O , is the escape time into the normal metal substrate. Secondly current flow also produces a drift of 

■ the order-parameter. Scaling equations for the superconducting gap and the current are derived 

' and are found to be consistent with previous phenomenological treatments as long as a temperature 

T ~ Teff is included. The current induced drift is found to produce additional corrections to the 
scaling which are smaller by a factor of 0{-^7^ — ), Ef being the Fermi energy. 
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QJ • I. INTRODUCTION 

U ■ 

I Quantum critical phenomena is the study of how a system loses long range order at T ^ as a parameter of 
^ , the Hamiltonian is changed^'^. The non-commutativity of position and momentum in quantum mechanics implies 
that the spatial and temporal fluctuations of the order parameter are coupled to each other at the zero temperature 
H , quantum critical point. The effect of temperature on a quantum critical point has some generic features^ such as, a 
I ' non-zero temperature produces dephasing or decoherence that cuts off divergences in correlation lengths and times, 
y ' Thermal decoherence also decouples spatial and temporal fluctuations causing a crossover from quantum to classical 
Q behavior. 

^ [ While quantum phase transitions for systems in equilibrium have been extensively studied, a much less understood 
' t issue is the effect of a nonequilibrium probe such as current flow on a system in the vicinity of a quantum critical point. 
Scaling theories exist which assume that the primary effect of a nonequilibrium probe is to produce decoherence or an 
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effective temperature. Thus nonequilibrium scaling relations are obtained by replacing temperature in the equilibrium 
^vq ■ scaling relations by the appropriate nonequilibrium energy scale^. A microscopic treatment to justify this and in the 
, process also identify the appropriate nonequilibrium energy scale is often challenging as this requires a treatment 
0^ ' that goes beyond a linear response Kubo formula calculation. Only a handful of such treatments exist for magnetic- 
\ paramagnetic —i^iii^, and superfluid-insulator/metal quantum critical point^ d°d^ . 

. In this paper we revisit the problem of non-linear effects, in particular the effect of a uniform current flow on a system 
' near a superconductor- metal quantum critical point. Existing studies have so far involved writing phenomenological 
OO , effective theories for a charged order-parameter in the presence of an electric field and/or external dissipation^ii^. 

' In this paper we carry out a fully microscopic derivation of the appropriate nonequilibrium effective theory starting 
^ , from a fermionic model under external drive. In doing so, we address the issue of how the fermionic system reaches a 
nonequilibrium steady state, and find that the underlying nonequilibrium fermions give rise to additional terms in the 
. _ effective theory for the charged order-parameter that were previously missed. We then proceed to determine the effect 
. of these terms on scaling near the equilibrium quantum critical point. (Note that the observation that nonequilibrium 
■ ■ ' electrons can significantly modify the scaling near critical points was also pointed out in^, and was experimentally 
observed in thin films of Bi ini^) . 

The geometry that will be studied (shown schematically in Fig [T]) is a 2d itinerant electron system where the 
electrons interact with each other via an attractive interaction. This system is driven out of equilibrium by an in- 
plane electric field so that a current flows through the bulk of the system. In addition, the 2d layer is coupled to an 
underlying normal metal substrate with which it can exchange particle as well as energy which thus serves as a heat 
sink that allows the layer to reach a nonequilibrium steady state. The coupling to the substrate also provides a source 
of dissipation that when made sufficiently large can destroy superconductivity in the layer—. (Note that the model 
in Fig [1] for the case of repulsive interactions of the electrons in the layer and for parameters that are such that the 
system is near a ferromagnetic-paramagnetic quantum critical point was studied in*). 

As we shall show, in equilibrium the effective action for the superconducting fiuctuations on the ordered and the 
disordered side has a local Caldiera-Leggett dissipation typical of systems where particle number is not conservedi^. 
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Effective theories for superconducting fluctuations witfi local dissipation have been extensively studied in equilibrium, 
14,15^16,17 ^ but hardly at all out of equilibrium (with the exception o^) . Out of equilibrium, our microscopic treatment 
reveals three effects of current flow, one is the usual minimal coupling of the current to the charged order-parameter, 
second is a source of noise which for low frequencies and long wavelength fluctuations of the order-parameter essentially 
acts as an effective temperature which equals the typical energy an electron gains on being accelerated by the external 
electric field. In our model Tefj = eEvpTsc, where e and vp are the charge and Fermi velocity of the electrons, and 
Tsc is the inelastic scattering time or escape time into the normal metal reservoirs. Thirdly, we also find that current 
fiow can cause the order-parameter to drift with a drift velocity vd = ^Tsc, m being the mass of the electrons. 

We briefly mention the relation of the work presented here to that 0^= which also had an extrinsic dissipation which 
was introduced phenomenologically. Thus in their model, current affects the order-parameter only via its minimal 
coupling to it, and the properties of the dissipative reservoir were unaffected by the current flow. In our model, the 
dissipation originates via the coupling of the superconducting order-parameter to the underlying normal electrons, 
whose properties are itself modified due to an external drive. Taking this effect into account shows that the order- 
parameter is subjected to a noise and also drifts with the current. As we shall show, in the quantum-critical regime, 
current noise gives corrections to the scaling which are of 0{{Tef fTsc)"^^"^), where T^ffTsc is the ratio of the typical 
energy gained from the electric field between collisions, and the energy lost due to inelastic scattering. Since, in 
our model Tsc is largely independent of the electric field (it may acquire some corrections at large electric fields), 
the corrections to scaling due to noise {{Tef jTsc)"^^^ <C 1) is subdominant in the quantum critical regime, with the 
dominant scaling behavior being that derived in^. In the quantum disordered regime however, direct coupling and 
noise effects are found to be equally important. Current drift on the other hand gives a correction which is smaller 
by an additional factor of 0{1/ EpTsc) where Ep is the fermi energy. 

The paper is organized as follows. The model is presented in Section HIl and is treated within a Keldysh path integral 
approach which will allow us to study out-of-equilibrium effects. We first study the equilibrium properties of the system 
by performing a mean-field treatment in section Hill which reveals a dissipation induced quantum critical point, which 
can also be understood as a proximity effect. A derivation of the effective action for the superconducting fluctuations 
about the equilibrium ordered side is presented in Appendix |^ and the origin of a local Caldiera-Leggett dissipation 
arising due to nonconserved particle number is highlighted. Fluctuation about the nonequilibrium disordered state is 
studied in section |V] and the new terms in the bosonic theory corresponding to current noise and drift are derived. 
Scaling equations for the gap and the current are derived in section I VI I Many of the details of the derivation have 
been relegated to the appendices. Finally we conclude in section IVIII where we discuss our results in the context of 
existing experiments. 



II. MODEL 



We consider a model of electrons in a 2d layer that interact via a short ranged attractive interaction responsible for a 
superconducting instability and are coupled via tunneling to a reservoir of non-interacting electrons. The Hamiltonian 
for the system is 

H = Hbath + Hiayer + Hiayer-bath (1) 

where Hiayer is the interacting electron layer whose critical properties we are interested in, Hhath represents the 
reservoir, while Hiayer-bath represents the coupling between the two. 

Hiayer = E (j ~ AV^j^t^X^T (2) 

Hbath = E '^L,fc,o-4,,fc,o-Cfc.,fc,(T (3) 



k^.k 

Hiayer — bath — 



a,kz^k 

a is the spin label, c represent the reservoir electrons, fc^ is the momentum transverse to the superconductor-bath 
interface and is not conserved on tunneling, while fc is the momentum within the layer. We assume the superconductor- 
bath interface to be smooth, so that the in-plane momentum is conserved on tunneling. The schematic of the model is 
shown in Fig[T] In addition the electrons in the interacting layer are subjected to a dc electric field which we represent 
via a vector potential A = —cEt. (Wc will set h = 1). 
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FIG. 1: A 2d itinerant electron system near a superconducting instability and driven out of equilibrium by application of an 
in-plane electric field. A steady state is reached via coupling to a normal metal substrate. 



We write the Keldysh action for this mode l^^'^^i^° , 



(5) 



where ± labels the Keldysh time-ordering, L'^'^^ is the action for the reservoir electrons, while L'^ is the Keldysh action 
for the layer electrons and the coupling with the reservoir. 



a 



dt 2m \ i c , 



C±a 



(6) 
(7) 



where /i is a chemical potential. 

We perform a Hubbard Stratonovich decoupling of the attractive interaction 

expUA / dtd''a;[*_T^'-i«'-i^'-T-*+T*+i*+i*+T] ) - / I?[A±,Ayexp 



exp (^i J dtdf^x (A_*_|*_x + A*_*_^*_| - A+*+t^+x - A;^'+x*+t; 



(8) 



and in the process introduce the bosonic fields A± which represent superconducting fluctuations. Using Nambu 
notation — ^^^^^ j ^± — (^±T V'±i) ,c± — {^^^ ~ (^±T c±j,) the Lagrangian becomes, 



~t 



V 

1 

-1 



9 _ 1 fY_ 

ot 2m \ I ^ J 



Ah 



c± + h.c. 



A 



(9) 



The electronic degrees of freedom may now be formally integrated out, resulting in a Keldysh action entirely in 
terms of the fluctuating fields A±. A rotation to retarded, advanced, Keldysh space leads toi^ 



Zk^ j V [A,,c;, a;,;] exp {TrXng-^) exp J dtd'' 



2A;A,, + 2A:,A, 



A 



(10) 



— — 2"^, Ac; — 1^^'^ are respectively the quantum and classical components of the fluctuating fields. 



where A^ 

5 is a 4 X 4 matrix in Nambu and Keldysh {r^.y^z) space which obeys the Dyson equation 

= Go' H 















) 









(11) 



where Gq is the exact Green's function for non-interacting electrons coupled to reservoirs and subjected to an external 

fG" G'^\ 

electric field. The full Green's function G may be written as follows in Nambu and Keldysh space, G ^ I q qa j 
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where the Q^'^'^ are the following 2x2 matrices, 



1 _ 1 /r^R pR 

and the retarded Green's functions arc defined as, 

G^'{x, t; x', t') = ~ie{t - t'){{^j^{x, t), (a;', t')}) = Gf (x, t; x' , t') (12) 

G«(a;, t; x' , t') = - t), V'i (a^', t')}) = -Gf (xY; x, (13) 

F'^ix, t- x', t') = -z9{t - t'){{i:^ix, t), (x', t')}) (14) 

F«(x, i; x', t') = -z0(t - t'){{4,iix, t), (a:', t')}) (15) 

and the Keldysh Green's functions are 

G^(x, i; x', i') = [V'tCa;, i), V^tC^^'. ^0] > = Gf{x, t; x', t') (16) 

G^{x,t;x',t') = -i([V;j^(x,t),i/'|(x',i')]) = -Gf (xY;x,i) (17) 

F^(x,i;x',<') = -z([V'T(a;,i),V^i(a^',0]) (18) 

F^(x,i;x',t') =-*([V^i(a:,i),V^T(2^''0]) (19) 

III. MEAN FIELD TREATMENT IN EQUILIBRIUM {E = 0) 

The mean-field equations may be obtained by minimizing Eg. 1101 with respect to the quantum (A^) and classical 
(Ac/) fluctuations of the order-parameter. A Ginzburg-Landau action is then obtained by expanding the Keldysh 
functional in fluctuations about the mean field solution. We first outline these steps for the equilibrium case i.e., 
A = Et = 0, before turning to the nonequilibrium case. 

In equilibrium the single particle Green's function Qq may be easily obtained. In Fourier space the retarded Green's 
functions are, 

G-^, = ^ - a - (20) 



'07?, 
"OR 

where = ^ fJ-, and the self-energies E arise due to coupling to reservoirs and have the form. 



Go-^ = ^+a--s« (21) 



- y — ? ^ -«r (22) 



Lu — — el + fi + iS 



^""^y -, — -, ? - -»r (23) 

where F = irpt'^ , p being the density of states of the reservoirs. We have taken the reservoir dispersion in Eq. Oto be 
^k, fc ^ ^fe, + ^k- Note that we will interchangeably use the notation 

Tsc = ^ (24) 

to represent the typical escape time into the reservoirs. 

For a reservoir in equilibrium at temperature T, the Keldysh self energies of the layer electrons due to coupling to 
the reservoir obey the fluctuation-dissipation theorem, 

S^ = -2irtanh— (25) 
2T ^ ' 

E^ = -2iFtanh^ (26) 
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Moreover, 







It therefore follows that 







iV 



=> Gro 



1 



{L0 + ^Tr-a 



-a 




iT 



and the fluctuation-dissipation theorem is obeyed so that, 



:IOK 



-2iT 



tanh ■ 



2T 



GorGoa 



(G 



OR 



:/0A 



tanh ■ 



2T 



(27) 
(28) 



(29) 



(30) 



The single particle Green's functions computed above for electrons coupled to reservoirs seem identical to those for 
electrons scattering elastically off impurities. However, the difference between the two systems will be apparent in 
the single particle level in the next section when an electric field is applied. In that case, while there is no steady 
state for electrons scattering off static impurities, the coupling to a reservoir in our model will be shown to provide an 
inelastic mechanism which will allow the system to reach a nonequilibrium steady state. The difference between the 
two systems is also apparent in equilibrium when electronic response and correlation functions are computed. For a 
disordered system appropriate disorder averaging gives answers which are consistent with a closed system characterized 
by conserved particle number. For our system the response and correlation functions (computed in Section |V} will 
refiect the fact that the system is open since electrons can escape into the reservoir. 



We now expand the Trln in Eq [TO] about Aq 



Ag + Aq, Ac/ 



Ac/, and determine Ag,Ao so that the 



resultant action is a minimum with respect to both quantum and classical fluctuations of the order-parameter. By 
choosing A^ = 0, the resultant Keldysh action is automatically minimized with respect to classical fluctuations of the 
order-parameter^^, whereas Aq will be determined by minimizing with respect to the quantum fluctuations. Thus the 
mean flcld Green's function is the matrix 



^mf ^ Go 



Ao 
Ao 



«) To 



where the retarded component is 



GRmf — 



and the Keldysh component is. 



GRmf 

1 



oj - £,k+iT Ao 

Ao uj + ^k + iT 

-Ao 



=> 



+ £.k + «r 



(w + tvf - Ag - f 2 V -Ao - 6 + «r 



[tanh ^] 



Glimf = — 2ir [tanh ^] GRmfGAmf ~ tauh ^ [GRmf ~ GAmf 
u)+^k+iT uj+^k—iF — Ao 



1 



-A 



Aq is obtained by minimizing Eq. [TO] with respect to A* which leads to the self-consistent gap equation, 

-2iAo 



A 



■ Tr [F, 







The above equation, together with Eg . implies the following condition for a superconducting instability 



7 > / — tanh — 



(31) 

(32) 
(33) 

(34) 
(35) 

(36) 
(37) 
(38) 
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For r = 0, r 7^ 0, the above is the usual equation for the BCS mean-field transition temperature. On the other hand, 
at T = 0, r 7^ 0, Eq. [35] yields a critical value of the dissipation r^, above which the system becomes normal 



1 1 ^bcs 

— r = In 



(39) 



where v is the single particle density of states at the Fermi energy. The destruction of superconductivity due to large 
r can be understood in a simple way as a proximity effect. F measures the amount by which the states in the 2d 
layer broaden due to hybridization with the normal metal leads. Thus the larger is F the more the states in the 
layer acquire the property of the underlying normal metal. For an early reference on similar proximity effect induced 
destruction of superconductivity in a thin superconducting film deposited on a normal metal see^^. 

A second interesting feature of this model is a gapless spectrum even when there are nonzero superconducting corre- 
lations. To see this we evaluate the quasi-particle density of states per spin direction N[ijj) = ^ {G^ — G^) (fc, uo), 
and find that at it is nonvanishing for lo < Aq. In particular at zero frequency 



N{lu = 0) 



(40) 



This appearance of gapless superconductivity is also a typical property of superconductor-metal interfaces^. 

To fully understand the equilibrium properties a derivation of the superconducting action on the ordered side is 
presented in Appendix \X[ To keep the notation simple, this derivation has been done for the partition function, and 
the effective action for the phase fluctuations in imaginary time is found to be (Eq. lAlip 



dr 



ci [droy 



where 



9_ 

27r 



5s 



(41) 



(42) 



(Eq. IA17[) . The distinguishing feature is the dissipative term Sy, which arises because the superconducting layer is 
characterized by non-conserved particle number. Since the physical quantity is the voltage fluctuations V = drO in 
the superconducting layer relative to the normal substrate, it is instructive to rewrite Eq. |42] after an integration by 
parts. 



Sy = -^J d'^x j dr j dT'dre{x,T)\n\T ~T'\dr'9{x,T') 



(43) 



IV. MEAN FIELD TREATMENT OUT OF EQUILIBRIUM 

We now turn to a mean-field treatment of the out of equilibrium current carrying case. A mean-field approach relies 
on the assumption that even in the presence of current flow the transition from the normal to the superconducting 
side is second order. In what follows the mean-field phase boundary in the current vs. equilibrium superconducting 
gap plane will be derived coming in from the disordered side. This will be followed by a discussion of the validity of 
me an- field. 

The self-consistent mean-field equations for the current carrying case is still given by Eq. I31[ however we now have 
to evaluate Qq which is the Green's functions for the layer electrons coupled to external reservoirs and subjected to 

an electric field^. An analytic solution may be obtained in the limit where l/EpTgc ^ 1 and T^ffTsc ^ 1, 

where Tsc = 1/(2F) is the typical escape time into the reservoir, and T^ff — eEvpTsc is an effective temperature that 
characterizes the steady state distribution function of the layer electrons. The details are presented in Appendix IB] 
Within this approximation the mean-field green's functions G'^^ remain unchanged and are still given by Eq. 1331 

while Qmf changes due to a nonequilibrium electronic distribution function. In what follows, both in the mean-field 
treatment of this section, and the study of fluctuations in the next section and in Appendix [Cl we will make the 
additional assumption that T^ffTsc ^ 1- 

Using the expression for the distribution function in Eq. IB15I IB16I IB17[ the self-consistent gap equation Eq [571 
becomes, 

4=-rf r^^fi-e-w^) fdc 



1 



[U! ■ 



iF)2 - ^2 



A2 
^0 



{lo - iF)2 - ^2 



A2 



(44) 
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We may approximate /^^ ^ J^'" ^ (^1 - e'V?^^ j ~ J^^'^'=^= Defining T^^e// as the critical value of the current 
induced effective temperature for which Ag = in Eq. I44[ and relating ^ to the gap A^q in equilibrium we find 



arctan — = In — bcs __ ^^t-^ 

An approximate solution to the above equation is Tc^eff — A^q/V^. Note that this current induced loss of order is 
a heating effect arising due to a highly broadened electron distribution function, and is not the same as the Landau- 
criterion for the critical current for breakdown of superfluidity. 

We now turn to the discussion of the validity of the above mean-field treatment. Firstly, within mean-field T^ff 
arises solely out of noise due to normal electron current, and leaves out the fact that current due to superconducting 
fluctuations also contribute to noise that can modify T^ff- Secondly, other scenarios for a current induced transition 
are possible. For example, on the superconducting side there is no dissipative current, so that Tg// ~ 0. In this case it 
is possible to have a supercurrent induced first order transition from a superconducting state to a normal or resistive 
state as discussed in^^i^i. Whether the actual transition is a second order heating effect as predicted by the mean-field 
treatment on the disordered side or a first order transition depends on whether the critical current for the first order 
transition is larger or smaller than the current corresponding to Tc.e//- In general this is a complex question that we 
do not address further. 

Instead, in what follows we will derive a nonequilibrium Ginzburg-Landau theory for the superconducting fluctu- 
ations on the normal side where most of the current is due to normal electrons so that there is a well defined T^ff. 
We will then use this to study how the gap and the current due to superconducting fluctuations scale due to electric 
field close to the quantum critical point and outside the fluctuation dominated Ginzburg regime. 



V. FLUCTUATION ABOUT NONEQUILIBRIUM DISORDERED STATE: DERIVATION OF THE 

KELDYSH GINZBURG-LANDAU FUNCTIONAL 

We now turn to the discussion of fluctuations about the mean-field disordered state ( Aq = in Eq. [nj) . In doing so 
we will also highlight the difference between how the electric field affects magnetic-fluctuations^ and superconducting 
fluctuations. 

Expanding the Tr In in Eq.[TO]in the usual way TrlnQ-^ = Tr InG^^ + TrGoA - iTr^oA^oA + . . . where 



^ = [a*, ) + [A*q o) ® 

one obtains an effective action 

Zk= jv [A,,,,, A;^,,] e-^Sl-^Sl-^S% + ... (47) 



with 



Si = Tr {Al a:,) [ 2 "gg, a +^gg\ ^ j (43) 



In position and time space 1 = x,t and 2 = x't' 

ngg(l, 2) = -^ [Gofl(l, 2)Gok(2, 1) + Gok(1, 2)Goa(2, 1)] (49) 
n^g(l, 2) = -^ [Gok{1. 2)Goif(2, 1) + Gofl,(l, 2)Goa(2, 1) + Goa(1, 2)Goi?(2, 1)] (50) 

Note that on the disordered side, terms cubic order in the superconducting fluctuations are absent (5*^ — 0), while 
Sjf has the form^. 

«.A;'Afp +C.C. (51) 



4=1. ..4 



We will treat only within a one-loop mean-fleld approximation. For this only the coupling constant ui will play 
a role. 
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Since Gfl,(l, 2) = -Ga(2, 1) Gnik, uj) = -G^C-fc, -w), Gk(1, 2) = -Gk(2, 1) Gx(fc, = -Gxi-k, -w), we 
may write 

Tr[A;nggAe,] = iTrA;(l) [Gofl(l, 2)Gox(l, 2) + Goa-(1, 2)Go/?(1, 2)] A,i{2) (52) 
Tr[A;n^gA,] = iTr A*{1) [Gok(1, 2)Gok(1, 2) + Gor{1, 2)Gob{1, 2) + Goa(1, 2)Goa(1, 2)] A,(2) (53) 

The coefficients 11 depend explicitly on the electric field. Since we use the Gauge A — —cEt, it is most convenient 
to go into momentum and time space, 

Tr[A;nggA,,] = iTrA*g{~q,ti) [GoBip + q;h,t2)GoK{-p]h,t2) + Gok{p + q;h,t2)GoR{-p;h,t2)] A,i{q,t^ (54) 

Note that for magnetic fluctuations, the above expression would have had the form^ 

Tr[mqIl^Qmci] = iTrm*{-q,ti) [Gor.{p + q;ti,t2)GoK{p;t2,ti) + Gqk{p + q;ti,t2)GoA{p;t2,ti)] mci{q,t2) (55) 

As shown in Appendix [Bl if the single particle Green's functions are written in terms of the canonical momentum 



k = p + eET where T = 
may write Eq. [54] as 



they become time translationally invariant. Thus for superconducting fluctuations one 



Tr[A;n«gA,,] = zTrA;(-q, ^l){GoJ^(p + q + eET- h - t2)GoKhp + eET; h - t2) 
+GoKip + q + eET; h - t2)GaR{~p + eET; h - t2)}A,i{q, is) 

In terms of the canonical momentum k = p + eET, the above becomes 

Tr[A;U^QA,i]=iTrA;{-q,h){GoRik + q;h-t2)GoK{~k + 2eET;h~t2) 
+GoK{k + q;h-t2)GoR{-k + 2eET;ti-t2)}A,i{q,t2) 

or shifting variables k k + 2eET one may write 

Tr[A;ng^Ac,] = iTrAli-q, ti){Gm{k + q + 2eET; h - t2)GoK{-k; h - ia) 
+GoK{k + q + 2eET;ti - t2)GoR{-k;h - i2)}A,,(g, ^2) 
Following the same steps for magnetic fluctuations we get 
Tr[mqIi^Qmci\ = iTrm*{-q,ti) 



(56) 



(57) 



(58) 



(59) 



Gor{p + q + eET; h - t2)GoK{p + eET, t2 - h) + Gok{p + q + eET; h - t2)GoAip + eET; t2 - h) mci{q, ^2) 



Rewriting the above in terms of the canonical momentum k = p + eET, one finds, 
Tr[mgU^(,m,i] = iTrm*g{~q,ti) 
Gofl(fc + q;ti- t2)GoK{k,t2 - h) + Gaxik + q;ti- t2)GoA{k;t2 - h) mci(q,t2) 



(60) 



Thus Eq. [SD] and [S5] highlight the difference between the coupling of the electric field to the magnetic and super- 
conducting order parameters. In Eq. [BDl all dependence of the electric field is via the modification of the Green's 
functions G^'^ at steady state, and there is no direct coupling between the electric field and the order-parameter. 
On the other hand Eq. 1581 depends on the combination {q 4- 2eET) which is the usual minimal coupling of the charged 
superconducting fluctuation and an external electric field. 

Thus to summarize, upto quadratic order, the Keldysh action for superconducting fluctuations in the presence of 
an electric field may be written as 



o2 



dti 



(A*(-<r,ti) A:,(-g,ii)) 



q + 2eET,ti-t2 



jSih -t2)+ q + 2eET, h - t2 



(61) 



}6{h - t2) + ( g + 2eET, h - t2 







'A,((f,t2)' 

Aci{q,t2), 
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where T = We now discuss the coefhcients Yl^'^'^ and highhght the appearance of new current dependent 

terms that were missed in previous phenomenological treatments. 

We expand the 11 bubbles in powers of (^q + 2eET^ to obtain, 



U-lih -t2) + E- {q + 2eET) nf (ti - ta) + {q + 2eET) U^ih ~ ta) 



(^q + 2eET;ti-t2^ = 

+ 2eET; h - ^2) - TIq (h - ^2) + C((g + 2eETf) 



It is convenient to Fourier transform the above expressions so that 



n'' q + 2eET: <i - i2 = 



^-in(ti-t2) 



2tt 



fioi^) + E ■ {q + 2eET^ Iff + {q + 2eET^ ^ 



(62) 
(63) 

(64) 



Each of can be evaluated as a power series in £7 (see Appendix [C] for details). Keeping terms to 

2 



0{Q.,\q + 2eETj ) one obtains, 

6ih - t2) + -^n^ [q + 2eET; h - = <^(ti - ^2) 
where, as derived in Appendix [Cl 

V = 

6 = 

7 = 




eE- (q + 2eET 



+ S + j(q + 2eET] 



i>Xt. 
1 

\v- 



X 



n^(o,o) 



65) 

(66) 
(67) 
(68) 



The first term on the r.h.s of Eq.[65]is the overdamped dynamics associated with non-conservation of particle number, 
while the second term is of the form V]j ■ ^9 + 2eAtj and represents current induced drift at velocity 

T.r_eE 



VD 



(69) 



The difference with^ is the appearance of the above drift term, along with a change in the noise properties of the 
reservoir (represented by H^) due to current flow. In particular, we find the following electric-field dependence of 11^ 
in 2d (see Appendix [Cl for details) 



^\+T,ff I ^Icos 



27r' 



(70) 



as opposed to a current independent cx in the model studied in^. 

Note that for a Id system, the structure of n^'"* remain the same as in 2d, while 11^^ acquires the form in Eq. IC6I 
Qualitatively it has the same structure as Eq. [TDlin that 11^ cx |0| when > T^ff, and 11^ cx T^ff when il = 0. 

We now turn to the evaluation of the gap-equation and the current due to superconducting fluctuations to see what 
role these new terms due to current induced noise and drift play. 



VI. EVALUATION OF SELF-CONSISTENT GAP AND CURRENT DUE TO SUPERCONDUCTING 

FLUCTUATIONS 

In order to derive the self-consistent gap equation and the fluctuation conductivity, as inS. we will work to quadratic 
order (Eq. I6ip . treating the quartic term in superconducting fluctuations (Eq. I5ip within a one-loop mean-fleld 
approximation. 

We may define the retarded, advanced and Keldysh component of the Green's functions for the superconducting 
fluctuations as follows 

Z?«(l,2) = -z0(ii -t2)([A(l), A*(2)]) = -z(A,(1)A:,(2)) (71) 
D^{1, 2) = ie{t2 - ii)([A(l), A*(2)]) = -z(A,Kl)A*(2)) (72) 
D^{1, 2) = -z({A(l), A*(2)}) = -*(A,,(1)A:,(2)) (73) 
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From Eqns [611 [65] and [TQl the equation of motion obeyed by the above Green's functions are 



where 



2eEti 



The above equation corresponds to overdamped dynamics and may be solved easily, 



where 



D^iq;ti,t2) ^ -9{ti-t2)-exp(-- f dr Egir) - iijVD ■ {q + ^eEr) 

V \ V Jt2 '- 



1 



D^{q;tut2) = -0{t2-ti)-ex^p 

V 



VJt, 



dr 



eg(r) + irjVD ■ [q + 2ei?rj 



eqir) = 5 + j (q + 2eETj 



(74) 
(75) 

(76) 
(77) 

(78) 



The self-consistent gap equation is 



A. Self-consistent gap equation 



S = <5o +Ml(|Ad|2) 



where ui ^ -^^^'^t'sc ^^^^ 

{\Aci\'^)=iD'^{x,t;x,t):^-i J d2d3D^{l,2)U'^{2,3)D^{3,l) 
d^q ' 



{2n) 



dh I dt2D'\q-tM)^'^{tl,t2)D''{q-M,t) 



In Fourier space n^(ti,t2) = / ^e-*'^(*i-*2)n^(f^) which together with Eqns [7S1I771 1701 give 



cl\ ) = 



d^q f dn 



(79) 



(80) 



{2n 



,2 / i^n^(") / dt. 



J_g-jf2(ti-t2)g--^ /4 dri[S+y(q+2eETif-invoiq+2eETi)] 



g~7;ft\ dT2[S+',{q+2eET2f+ivVD-{q+2eEr2)] ^^^^ 

Changing variables to the canonical momentum k = q + 2eEt, the explicit dependence on t goes away, and one obtains 



^tn(x-y)+ivD-k{x-y)-iVD-eE{x^-y^) ^~^{5+'yk^){x+y) ^-^{2eEf{x^+v'>) ^-^k-^^^ 



(82) 



It is convenient to perform the momentum integrals, which gives. 



(|Ac 



r Jo 



dx / dye 



iVt{x ~y) —iv u -eE^x^ ~y^) 



1 



x + y 



-^ieEfix'+y')JieEf^^^^{l+'^^(^}^ 



e 



(83) 



After this the manipulations are similar to^. It is convenient to change variables to u = x + y, v = x — y, so that 
/o°° /o°° = 5 /o°° ^" I-u giving 



47r2 7 7 27r ^ ' 2r)^ Jq u 



dv 



.^Sn-eEuv -^u -^{eEruiu^+3v^} i^ieEf^^^^{l+^^^y 



(84) 
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Now we approximate the expression for II-^ in Eq. [TOlas {Q) ~ 2i?] \Q\ + ^^^^9{\n\ — T^ff) 
dimensionless variables Cl ~ Qt], u/i] ^ u,v/ri v 



We also define 



Tejf = TeffT] 

E = eErjy/j 



(85) 
(86) 



in terms of which 



(|Ac/P> 



1 



- IT tf 

\n\ + ^lle{\n\-T,fj) 

TT 



e2 („2^^2)2 



— / dvcoslilv E uv -\ E - 

u Jq \ 7m 27m u 



(87) 



Changing variables to w v/u, one gets 



(lAc/l^) = 



1 



- 2T 

|f^| + ^0(|l]|-re//) 



du I dv COS ( nuv - —E^u^ + ^^E^v 



7m 



27m 



g-5«g-^(l+6-'-3-^)-4i^™'(l+-') 



(88) 



It is now straightforward to see the role played by current drift. This term always arises in the combination ^E'^. 
Using Eq. [68l one finds it to be O ■ we shall show, the electric field scaling due to the current noise term is 

0{E^) in the quantum-disordered regime. Thus the drift gives corrections to this result by an amount which is smaller 
by a factor of ^ <C 1 {Ep = fi). Therefore in what follows we will drop the drift term from further analysis. 
Substituting Eq. 1881 in Eq. [751 adding subtracting terms one gets the following self-consistent gap equation 



27r7?7 



Ui 



dn. 



— ]n\ I du I dvcosinuvje 
27^777 J_oo 27r 



— — |f2| / du I dvcoa{nuv)e 



-5u 



12 (1 + 6^ ) 



dn 



27r7?7 J_oo 27r 



2T, 



e 12 -t- e 12—1 
du I ' dvcos (f]uv)e-*"e-^(i+^"'-3"') 



(89) 



(90) 



We introduce a frequency cutoff A in the first term in the above equation, and perform the frequency integral to 
obtain, 



Ul , 1 
■ In 



27r^777 6 



5o 



27r7?7 

Ml 



dn 



1 



— \n\ I du I dvcos {nuv) e^^'' 







tl+6V-3v^ 



_ £2 „3 ^ g2^3 

e 12 -I- e 12—1 



dn 



2vr777 2tt 



2T, 



eff 



0{\n\ - T^ff) 



pOO nl _2 3 

/ du dvcos {nuv) e~'^e~^i'+'^^"-^^") 
Jo Jo 



(91) 



(92) 



The remaining frequency integrals above are performed by introducing a cutoff e ^ l'^' in the argument, and then 

(A^2_|_„2j,2p U^ll^ 



setting A ^ = 0. For example, one integral evaluates to Lt/^-i^Q ricos(r2uu)e ^^/^ = — 7:^, while 



another is Li^-i^o /o°° ^in (f7u)e ^/"^ = i. 

In addition, by defining, 6^ — nj + [Aj — (5 In A] as the renormalized distance from the QCP, and by using 

27r2-,, 

In 1/(5 ^ 1, the self-consistent gap equation becomes. 



(51ni = ,5^- /" du \ dv- 



1 



%^(l+6"^-3^^)-p-4# 



(93) 



-I- / du—e 
u^ 



e2,.3 

e 12—1 



+ (^) du dv^sin{T,ffUv)e-'-e-^('+'^'-'^') 



1 

uv 
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The first three terms on the r.h.s was derived in§., whereas the last term is new and arises due to current noise and 
reflects the modification of the underlying electron distribution function. 

We discuss the solution of the gap equation in two regimes 
A. Quantum Disordered Regime, 5^ ^ E. Here Eq. [^Slcan be perturbatively expanded in powers of E to give 



E^ 



^eff 



h^+3(5«/ln3^)2 ' TT (5^7 In 



(94) 



While the first two terms in Eq. [94] were derived in^^ the last term is the correction due to current noise which 
essentially acts as an effective temperature. As discussed before, current drift will correct this result by a factor of 
0{1/EfTsc). 

B. Quantum Critical Regime, 6^ ^ E. Here one may set e"^" = 1 in Eq. [93l which in terms of a rescaled variable 
u — uE^^^, may be written as 



Sin- 



1 



^2/3 



/ du^ (e 
Jo V 



1 



du i dv- 



_ =1(1+61,2 



/ 2T' 



eff 



\^7r^2/3 

Defining the following functions 



du 



y = 



1 



31/324/3 



r ; 



dv- 

y 



(l + 6i;2-3i;4)'/'-(l + «2) 



du 



0.1165 



1.603 



(95) 



(96) 
(97) 



we find, 



{2Ef'y 

In 1 



(98) 



Again the first term was derived in^, while the second term above is the correction arising due to the modification 
of the distribution function of the underlying electrons. Using the definitions Eq. [551 [Ml and the expressions for 

^2 2/3 

and 7 in Eqns.[ni[nSl ^ - {TeffTsc) ' < 1. Thus this term only gives rise to subleading corrections within the 
model presented here where r^c is independent of the electric-field. In the conclusions we discuss the case of systems 
where t^c may have a strong electric-field dependence, and can in particular diverge as ^ 0. In this case it may be 
possible for the second term to dominate over the first. 



B. Expression for the current due to superconducting fluctuations 



We now turn to the evaluation of current due to superconducting fluctuations. The expression for the current is 
given by 

Changing variables to the canonical momentum k = q + 2eEt, and using Eq. 1741 we obtain the expression 

J = -1^1 J J ^U^{n)^ J^"^ dx J^^^ dye^''^'^-y^e-H'+'<''>^^^^ (100) 

Performing the momentum integral, one gets 

47r2 7 \ y J 2tt If Jq Jq [x + yf 
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As before we perform a change of variables to previously defined dimensionless variables to obtain, 



2e^E_ 
~2tt 



2^ 



11^1 



2T, 



elf , 



M-Teff) 



du I dvcos(n™)^(l + w2)e-'5"e-^(^+^'''-3"') (102) 
2 



Adding and subtracting terms in Eg 11021 



2e' E 



dn. 



+ 



h 27r 



2e2 E f2T, 



27r' 



du 



dv cos (riuu) — (1 + w^)e 



(103) 



h 27r \ TT 



dn 

2^ 



[^(1^1 - T^ff)] / du 



dv cos {nuv)^{l + „2)g-5ng- ^^^(1+6.^-3.^ 



One of the integrals that may be performed is Lt 



dnne--"/^-^ [uficosiun) + (-1 + u^n^)sm{un) 



Lt 



0°° dflile-^^^ /o^ cos(l]uw) (1 + ^2) 

= 0. Moreover using 



7t/2 
A 



LiA-i_o /o°° COS (f]ut;)e~^V'^ 



e^E r°° 



du 



e^E f2T, 



1 



^ the expression for the current becomes 

g-4#(lW-3.-)_g-4# 



dv ^ ' 



(l + i^^)e- 



(104) 



1 - M 

/ dv— sin (Teffuv) -{I + v')e-''^e 







As before we discuss the following two cases 
A. Quantum Disordered Regime, ^ E. In this regime we find 



J = 



e'E 



4 f 2 



15 Stt (52 



(105) 



The second term above is the correction to the results o^ due to the effective temperature of the nonequilibrium 
electrons. 

B. Quantum Critical Regime, 5^ ^ E. Here we obtain the result 



J = 



e'E 



^^{l±v^ In (1 + 6^;2 - 3^;4) + 



du / dw(l + i;2)e-^(i+6-'-3.^) 



Computing the above integrals we find 



0.46e2£; 
h 



1 + 0.82 



f2 
_eff 

£;4/3 



(106) 



(107) 



The first term is the universal conductivity found in^, while the second term is the contribution due to current noise. 
As discussed after Eq. [§51 this correction is of O {{T^f fTscY^^) and is therefore subleading for this model of electric 
field independent t^c- 

It is instructive to see how the current due to superconducting fluctuations in the quantum critical regime get 
modified for a Id system. The steps in the derivation are the same except that there is only one momentum integral 
in Eq. [lOOl We find 



Jld — 



(108) 



where g{v) = 1 + 6v^ — 3v*. 

Eq. 11081 shows that unlike 2d, the response to the electric field in the quantum critical regime is highly nonlinear, 
with Jld cx E^/^. Current noise here too gives subleading corrections of 0{{Tef fTsc)'^^^). 

The results presented above are for the case of (5 > 0, i.e, the system is on the normal side in equilibrium. The case 
of 5 < and large electric fields so that one is on the current/supercurrent induced disordered side can be analyzed 
by employing a purely classical Ginzburg-Landau theory corresponding to a temperature T — T^ff. The computation 
of the non-linear response would foUow^^, where using their results one expects the fluctuation current in dimension 
d to be, Jd oc ^'[iJI)/3 E. 
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VII. CONCLUSIONS 



In summary starting from a fermionic model under external drive, we have presented a microscopic derivation of the 
effect of current flow on a superconducting order-parameter. Our microscopic treatment reveals that current besides 
directly coupling to the order-parameter also produces a noise and a drift of the order-parameter, the origin of which 
is the underlying nonequilibrium electron gas. We study the effects of these new terms on scaling near the equilibrium 
quantum critical point. Scaling equations when only the direct coupling between the order-parameter and the electric 
field is present was derived by Dalidovich and Phillips^ in a phenomenological approach. Here we find that current 
drift gives a small correction of 0(1/ EpTgc) to their result. Current noise on the other hand gives corrections that 
are of C'((re//Tsc)^^^) in the quantum critical regime. In our model where t^c is independent of the electric field, this 
correction is subdominant to the effect of the direct coupling between the order-parameter and the electric-field. In 
the quantum disordered regime however the noise and direct-coupling effects are found to be equally dominant. 

One may easily imagine a scenario where noise effects dominate over direct-coupling effects both in the quantum- 
critical and quantum disordered regime. This would occur when Tsc ~ ^^e// ^^'^^^ p > 1, a physical situation for this 
being when the dominant inelastic scattering mechanism is due to phonons. There are several experiments involving 
electric-field scaling in thin films near a superconducting transition 1^,26^27 ^ Jigcussed irJ^., the results of many of 
these experiments can be explained only when taking into account noise effects due to a nonequilibrum electron gas. 
For example p = 2 for electron-phonon coupled MoGc thin films, clearly making T^ffTsc ^ 1 in these systems. 

Our derivation is valid on the normal side and outside the Ginzburg regime. Extension of the results of this paper 
to the nonequilibrium ordered side is currently in progress. 
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APPENDIX A: EFFECTIVE EQUILIBRIUM ACTION FOR FLUCTUATIONS ABOUT THE ORDERED 

STATE 



In order to understand the fluctuational properties on the ordered side in the absence of an applied electric field, 
the action will be derived for a partition function. 



Z = V[A, A*] exp 



lAr 



Tr In Q 



(Al) 



where in terms of a complex A = Aqe'^*^, A* = Aqg 

/ 



~ 2m \ I r I ' 



Aoe 



2i9 



V 



(A2) 



E = E are the self-energies due to coupling to the underlying metallic substrate with 

TT \ T 



(A3) 



The action may be written as an expansion in fluctuations in the magnitude Ao and phase 6 of the order parameter. 
In what follows we will consider only fluctuations in the phase as the fluctuations in the magnitude of A are gapped 



in the ordered phase. To this end it is convenient to introduce the unitary matrix U 
the Green's function as 









V 



Ao 



and transform 

(A4) 
(A5) 
(A6) 
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where (j) = drO, A = A — -VO and we have spht the above terms as follows 



5o' 



ii^n - ik+ iTsgn{uJn) Aq 

Ao iujn + Cfe + iTsgn{uj„) 



Xia = -i<73(i 

Xib = S — e 
X2 = -fa: 



Expanding to quadratic order in the fluctuations, the action for the superconductor takes the form 

2" 



Z = / P6'exp 



drd'^r 



ci (dref + C2 (ve - -A 

\ c 



Si 



(A7) 

(A8) 
(A9) 
(AlO) 



(All) 



where the first two terms above are the usual ones that arise in any superconductor with the coefhcients changed due 
to coupling to an underlying substrate. In particular, 



Cl 



C2 



1 1 



Tig _ 1 1 

2^ ~ 2m?djL- 



-^Tr [p^Go'^oGo'^o] , d — dimension 



(A12) 
(A13) 



The new feature is Ss which arises specifically due to coupling to external normal metal reservoirs and reflects the 
lack of gauge invariance associated with the non-conservation of particle number in the superconducting layer. To 
leading order in the fluctuation of the phase, 



5s - Tr [GXi 



[GXu] - Tr [(G(xi, xt')E{t\ t) + E(i, t')G{xt' , xt)) (l - e-^«(-^*')e»»(-.*) j 



Evaluating the above trace we obtain a Caldiera-Leggett type local damping, 



ie{x) 



where 



3 = 4r 



,x/AfTT^, 



Fourier transforming Eq. IA15I one gets. 



^ ^ I d'^x I dr I dr' 



ix,T)~e{x,T') 



(AM) 



(A15) 



(A16) 



(A17) 



APPENDIX B: DERIVATION OF STEADY STATE SINGLE PARTICLE GREEN'S FUNCTIONS 

1. Derivation of retarded Green's functions 



To obtain the retarded Green's function in the presence of an electric field and coupling to an external reservoir we 
need to solve the Dyson equation. 



[idt, - Ho{ti)] G^ih - t2) = 5{ti - h) + S^G^ 



where Hoit) = J2k^a ^ (^-i 



£4^)Vi V-E ,i*=-cMand 



E«(fc 



k^,k^ 



(Bl) 



(B2) 
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For energy independent tunneling amplitude, density of states, and using the fact that e\ j.^ = + the above 
expression simplifies to give an energy independent self-energy 

I]«(c.) = -27rfp J delSiu - el - ) = (B3) 

The above implies 

S^(<i,i2) = -*r<5(ii-i2) (B4) 



Substituting the above in Eq. IBli it is straightforward to show that the retarded Green's function in the presence of 
an electric field and coupling to leads is: 

where r = ti — ^2 and k = p + eET, (where we set h ~ 1). The above time integral in the argument may be performed 
to obtain the following series expansion 

G^(fc,r) = -i0(r)e-"^"-'i^(^^-A)'^^+-e-r" (B6) 

Now we define 

Teff = eEvpTsc (B7) 

r-i = 2r (B8) 



and Ep = vp/a with a being the lattice spacing. Then, the second term in the argument of the exponent in Eq. IB6I 

is {eEf{T^^d^ek/{dP)) = ^^"'sf^^ ^^g^r"^ < 1 and therefore may be neglected. A similar argument applies to the 
higher order terms. 

Thus, we may approximate the retarded Green's function by its value in the absence of an electric field, 

G^(fc, t) = -i9{T)e-""''' e~^'' (B9) 
provided k is chosen to be the canonical momentum. 

2. Derivation of steady state Keldysh Green's function 

The equation of motion obeyed by the Keldysh Green's function is 

'idt.-Ho) G^{hM) - 1 + S^G^' + S^G^ (BIO) 



ihM) {~id t, ~ Ha) =1 + G^E^+G^S^ (BU) 
Taking the difference between the equations IBIOI and IBlli one obtains 

{idt, + idt,)G^{hM) - e{h)G^{h,t2) + e{t2)G^{hM) - S^G^ + S^G^ - G^S^ - G^^S^ (B12) 
The solution for G^ may be obtained by using the ansatz 

Gn = G^ fK - JkGq (B13) 
where 1 — 2f — fx, with / the generalized distribution function. The equation of motion for is 
.dfK , .dfK , , 

^-Q^+^-Q^-^p—tA{ti)JK + ^p-^^A{t2)J -Jk + Jk-^ +^ =0 

Y^R _Y^A _ g^j-^j _ (S^ — I]^)(l — 25), g being the distribution function of the substrate. Fourier transforming 

Eq. IB14I with respect to the relative time t = ti — changing variables to the canonical momentum k ^ p + eET 
and expanding in E one finds that the distribution function at steady state obeys, 

. ^ + |/ f ,g . ^) + 1 f / f . i.) _ _L + ,1 (B„, 

dk duj\ dkj 2Aduj^\ dkj Tsc 
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The usual quasiclassical arguments imply that the first term in Eg IB14I is negligible while in the the weak field limit 
the third term may be dropped. With these simplifications wc find 



where 



/ = r + r 



(eE ■ VkTs, . 



2\l (^eE ■ f/fcTsc) 



(B15) 

(B16) 
(B17) 



where x = lo — ji and Vk = dek/dk. Substitution of Eas [BT6llB17l into Eq |B14l then shows that the neglect of the third 
term in Eq|B14jis justified when the coupling of the layer to the substrate is sufficiently weak ( EfJdk'-^ 
while the first term is negligible in the weak field limit T^ff ^ ^p/i Sdk'^ )• 



APPENDIX C: EVALUATION OF THE POLARIZATION BUBBLES WHEN E^O 



The retarded and Keldysh polarization bubbles may be expressed as an expansion in ( if + 2eET ) and ft as shown 
in Eq. EH In particular, 



U"{q + 2ET = Q,n) 



cLu 



J2 I ^ GoB.{k,uj + n)GoK{-k,-uj) + GoK{k,uj + n)GoR{-k,-uj) 



n" {q + 2ET = o,n) = 

^ \GoK{k, Lo + 17)GoA'(-fc, -c^) + GoR{k, oj + n)GoR{-k, -lj) + GoA{k, lo + n)GoAi-k, -^l 



(CI) 



(C2) 



Using Eq. Hi EHEq. EISlEq. iBTBllBTTl we find 



Im 



T 2^ 



{~2try 



' ' eE ■ vpTs, 



where 2r = t^J'. For eEvpTgf, < I/tjc, the above expression simplifies to 

Im [U^in)] = -2iVlvTsc 

In the same way, one finds 



^\+Teff I :r^|cos0|e 



2tt' 



(C3) 



(C4) 



(C5) 



It is also instructive to derive the expressions for the polarization bubbles in Id. While n^"^ have the same structure 
as in 2d, the noise 11^ has the form. 



\n\+Teffe -^ff 



(C6) 



Note that the above expressions were derived and used to study the effect of current flow on magnetic fluctuations 
in^. 
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In order to evaluate the coefficient of expansion in powers of q + 2eET in Eq. [65] we first write the expression for 
the polarization bubble in momentum-time space 

Il"{q + 2ET,ti-t2) = 

i J2 [GoR{k + q + 2eET, h - t2)GoKhk, h - ^2) + Gok(^ + q + 2eET, h - t2)GoR{-k, h - ^2)] (C7) 

k 

For convenience, we shift variables so that q + 2eET appears in the argument of the retarded functions so that 

i J2 [GoR{k + q + 2eET, h - t2)GoK{-k, h - ^2) + GoK{k, h - t2)Gofl(-fc + q + 2eET, h - is)] (C8) 



Using Eq. IB5[ the retarded Green's function can be expanded in a power series in Q = q + 2eET as follows, 



GoR(k + Q,T) = GoR{k,T) 



, k-Q .Q' I Q-k\ , 

1 — I r — I T T" 

m 2m 2 \ m I 



(C9) 



In the above we assume quadratic dispersion. Note that in equilibrium, the linear in Q term docs not survive the 
angle integration. On the other hand, a non-zero current picks a preferred direction so that this term for our case 
will no longer be zero. As we shall show, this term will give rise to current-drift. 
Fourier transforming Eq. IC9I with respect to r we get 



GoR{k + Q,uj) 



1 



k-Q d d 1 Q-k\ 



m doj 2m duo 2 \ m j duP' 
Thus Eq. ICIOI and IC8l lead to the following for the particle- hole symmetric case. 



(CIO) 



A(n^(g,o)-n>,o)) = ^g^-^ff^Jr. 



(CU) 



where 



— [ [ dijjsgn{uj) 

IT m J J 



{u - a)3 - 3r^(c^ - 6) 



(C12) 
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